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In nematic liquid crystals light is strongly scattered from director fluctu- 
ations. We are interested in the limit where the incoming light wave is 
scattered many times. Then, the light transport can be described by a dif- 
fusion equation for the energy density of light with diffusion constants D\\ 
and D^, respectively, parallel and perpendicular to the director. We start 
from a radiative transfer theory, connect the diffusion constants to the 
dynamic structure factor of director fluctuations, and shortly discuss our 
results. Temporal correlations of the diffusing light probe the dynamics 
of director modes on much shorter time scales than single light scattering 
experiments. To account for the decaying temporal correlations, one has 
to add an absorption term to the diffusion equation, which we also link to 
the dynamic structure factor. 
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1. INTRODUCTION 

More then one decade ago the discovery of coherent backscattering or weak 
localization of light in colloidal suspensions [1] initiated a tremendous re- 
search activity in multiple light scattering [2]. Researchers were attracted 
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by the possibility to achieve the equivalent to the Anderson localization of 
electrons in disordered solids [3] . So far strong localization of light waves 
has not been observed. 

In the theoretical description of coherent backscattering the diffusion 
limit for multiply scattered light is employed [4]. Photons are consid- 
ered as random walkers with a scattering mean free path /, measuring the 
length they travel between two scattering events, and a diffusion constant 
D — c/*/3, which involves the transport mean free path I* = — cosies). 
It stands for the path length beyond which the direction of propagation of 
a photon is randomized. The angular brackets denote an average over all 
possible scattering events and ■j^s is the scattering angle. In 1987 Wolf and 
Maret discovered that diffusing light could be used for spectroscopy [5], 
which was later called Diffusing Wave Spectroscopy (DWS) by Pine et al. 
[6] . This was an important step forward, since so far turbid systems could 
not be investigated with conventional dynamic light scattering. In DWS 
temporal correlations of the detected intensities decay much faster than in 
single scattering since phase shifts of electric field waves from many scat- 
tering events are added up. DWS, therefore, detects dynamic phenomena 
at much shorter time scales then normal dynamic light scattering [7]. 

Research on diffusing light and DWS has focused on isotropic me- 
dia, like colloidal suspensions [7, 8], emulsions [9], and foams [10]. Re- 
cently, diffusing photons were used for the imaging of objects [11, 12] 
which found its application in medicine [13]. In magnetic fields a pho- 
tonic Hall effect was discovered [14] . Interest in diffusing light in nematic 
liquid crystals was again initiated by the observation of coherent backscat- 
tering first in multidomain samples [15] and then in a perfectly aligned 
nematic state [16]. In treating multiply scattered light in the nematic 
phase, one has to deal with the anisotropy in light propagation and a 
different scattering mechanism. Light is not scattered from local objects 
like particles in colloidal suspensions but rather from long-range corre- 
lated director fluctuations [17, 18]. The theory for diffusing light in ne- 
matic liquid crystals was developed independently by Stark and Lubensky 
[19, 20, 21] and Tiggelen, Maynard, and Hciderich [22, 23]. For a review 
see ref. [24]. Its final content can be summarized in an anisotropic diffu- 
sion equation with absorption for the electric field autocorrelation function 
W{R, T, t) = {E(R, T + t/2) ■ eoE(R, T - t/2)): 
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§f-D\\V\-D^Vl-rl^{t) 



W{R,T,t) = 0, (1) 



where D\\ and D_l denote, respectively, the diffusion constants parallel and 
perpendicular to the nematic director. The dynamic absorption coefficient 
/i(t), measured in DWS experiments, governs the temporal decay of the 
autocorrelation function [fi{t = 0) = 0]. The mentioned theories link 
the parameters D\\, D±, and nit) to the structure factor of the director 
fluctuations. The absorption coefficient appears as an angular average 
over all director modes. Careful experiments by Kao et al. demonstrated 
the anisotropic diffusion and the application of DWS [21, 25] . Their results 
were in excellent agreement with theory. 

The developed theories used a Green function approach to arrive at the 
diffusion equation. In the present article we present a derivation without 
the extended formalism of Green functions but which nevertheless con- 
tains the main ideas. To do so we start from a radiative transfer theory 
which was first introduced as early as 1905 by Schuster to describe the 
transport of hght through the atmosphere [26]. The radiative transfer 
theory is basically equivalent to a Boltzmann equation for the energy den- 
sity of hght [27, 28]. Any coherent superposition of electric field waves 
are omitted. For anisotropic media it needs some refinements which we 
present in section 3. The advantage of the radiative transfer theory is that 
it only requires some intuitive understanding. It can, however, be derived 
directly from Maxwell's theory. For nematic liquid crystals this was done 
first by Romanov and Shalaginov [29]. 

The article intends to stimulate further research on diffusing hght in 
liquid crystals. There exists a wealth of material which scatters light very 
strongly, namely porous media filled with nematics [30] , polymer dispersed 
liqTiid crystals [31], focal conic textures in cholesterics [32], and the Blue 
Phase III. Also applications to lyotropic and polymeric liquid crystals and 
to liquid crystalline colloids could be of interest. 

Section 2 first summarizes light propagation in a homogeneous ne- 
matic phase emphasizing the differences to isotropic systems. Then we 
review single light scattering from director fluctuations and introduce the 
dynamic structure factor. Section 3 presents the transport equation for 
radiative transfer. An approximation for the diffusion constants is derived 
and the formula for the dynamic absorption coefficient is developed. The 



H. STARK 



main features of both quantities are discussed. 

2. LIGHT PROPAGATION AND SINGLE SCATTERING IN 
NEMATIC LIQUID CRYSTALS 

There are two necessary ingredients to treat multiple light scattering in a 
turbid medium: the propagation of a plane wave between scattering events 
in an effectively homogeneous system and the single scattering itself. We 
start with the first point. 

A. Light Propagation 

A homogeneous nematic with the equilibrium Prank director tiq represents 
a uniaxial and therefore birefringent material. It posseses the dielectric 
tensor [17, 18] 

£o = £±1 + Asrio <S>no , (2) 

where e\\ and e± denote the dielectric constants for electric fields, respec- 
tively, parallel and perpendicular to the director, and where Ae = e\\ — £_l 
stands for the dielectric anisotropy. There exist two characteristic light 
modes which we describe by plane waves for the electric field: E{r,t) = 
E°'ea{k) exp[i{—u!t + k ■ r)]. They are characterized by the polarization 
vector ea{k) and the refractive index na{k) = ck/u [33]. The ordinary 
light wave {a — 2) behaves as in an isotropic medium. It has a constant 
index of refraction n2 — ■•/sT and the polarization vector 62 (fc) is perpen- 
dicular to both Uq and k. For the extraordinary mode (a — 1), however, 
ni{k) depends on the angle -d between k and tiq: 

1 sin^f? cos^'j? , , 
^ = + . (3) 

The vector ei(fc) lies in the plane defined by tiq and k but is not perpen- 
dicular to the wave vector k. It makes sense to introduce the polarization 
vector d"{k) = eoea{k) for the electric displacement field. Then, the 
biorthogonality relation d"" ■ = 5^ holds [20, 34]. The normalization of 
the polarization vectors is chosen such that the energy density E -SQE/iTT 
of a light mode, averaged over on time period, becomes W"' = \E"'\'^/87r. 
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Besides the phase velocity Vpa = ck/ua there exists the group velocity 
Vga — Vfca;Q,(fc), where u!a{k) — ck/na{k) stands for the dispersion re- 
lation. Electromagnetic energy is transported along the direction of Vga- 
For ordinary light Vg2 — Vp2. For extraordinary light the group velocity 

/ cos-i? sin-i? \ 
Vgi = cni{k) no H ui (4) 

differs in direction and magnitude from the phase velocity. We introduced 
the unit vector iii perpendicular to both Uq and 62 (fc). In a medium with- 
out absorption the energy transport is described by the Poynting vector 
c{E X H)/A7i. It can be rewritten so that for each light mode it equals 
the energy density times the group velocity: = W^Vga [33]. 

B. Single Scattering 

Single light scattering from thermally activated director modes has been 
well understood for a long time [17, 18]. The fluctuating part of the direc- 
tor, 6n{r,t), induces fluctuations in the dielectric tensor, 5e = A6[6n ® 
Uq + Uq ® Sn], which scatter light. The typical scattering experiment in- 
volves incoming light with wave vector k°' — coUak / c and with polarization 
ea{k) which is partially scattered into light with wave vector — ujupq/c 
and polarization e0{q). In the weak-scattering approximation, the tempo- 
ral autocorrelation function of the scattered electric field is proportional 
to the dynamic structure factor 

Bap{k, q, t) = -^{Ssa^iqs, t) 0) ) , (5) 

where qs = — k"' denotes the scattering vector and S€a^{qs,t) = 
6/3(5) ■ ^^{Qs,t)ea{k) stands for the projection of 6e on the polarization 
vectors. In its final form we obtain for the structure factor: 



B^p{k,q,t)^{Ae)'—Y,N{a,P,S) 



4 2 



(5=1 



'Ksiqs) 



exp 



Ksjqs 



(6) 



The term k-oT / K^^qs) represents the thermally activated director modes. 
It involves the elastic coefficient Ks{qs) = K^q^ + Ka^y + A^i?^ with 
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Prank constants K^, magnetic field H, magnetic anisotropy Ax, and the 
components of Qs — {q±, q\\). The exponential factor reflects the diffusive 

natTirc of the director modes. The relaxation frequency depends on the 
viscosity rjs{qs) = 7 — fJ^{<ls), where the rotational viscosity 7 plays the 
important role. Finally, N{a, f3, 6) stands for a geometry factor which 
forbids ordinary-to-ordinary scattering and forward scattering along the 
director. From the structure factor we derive the scattering mean free 
path la{k) [20, 22, 35]: 

= E / -^^Bapik 0)4(g) , (7) 

which depends on the polarization a and the direction k of light. In the 
photon picture it gives the average path length a photon travels between 
two scattering events. For light intensity it gives the distance after which 
the initial intensity Iq has decayed to lo/e. However, multiple scattering 
events are totally neglected in the last interpretation. For a careful inter- 
pretation of the scattering mean free path in connection with k and Vga 
see ref. [21]. The mean free path is discussed in refs. [24, 35]. We only 
note, that /i(A;) goes to zero when H ^ since the structure factor for 
extraordinary-to-extraordinary scattering diverges for Qs 0. 

Finally, we make some comments about the symmetry of the structure 
factor. In isotropic systems B{k ■ q) only depends on the scattering angle 
•ds via k ■ q — cos'&g- Therefore, if expanded into spherical harmonics 
yi'm'{k) and Yi„i{q), it is fully diagonal: {l'm'\B{k ■ q)\lm) oc 5ii/5mm'. 
The structure factor for the director modes just posseses the rotational 
symmetry around the equilibrium director Uq. It can be written as a 
function of the relative azimuthal angle (p = ipq — (pk between q and 
k. Using the additional symmetry Bai3{k, q,0) = Bai3{—k,—q,0) the 
following expansion holds: 

B^^ik ^, 0) = E Kfsi^k: ^) cos[m{<p, - <pk)] , (8) 

m>0 

where '&k and '&q are polar angles of k and q with respect to riQ. It 
expresses the fact that {l'm'\Bai3{k, q,0)\lm) oc 6jnm' is diagonal in the 
index m but not in I. 
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3. RADIATIVE TRANSFER THEORY AND DIFFUSION AP- 
PROXIMATION 

In the following we deal with the temporal auto correlation function of 
the electric field 

nl{k)W^{R,T,t) = {E'?{R,T -t/2) EI*iR,T + t/2)) . (9) 

For t = 0, it stands for the energy density of a light wave at time T 
and space point R travelling into direction k with polarization a. The 
light frequency u is omitted. We pulled out a factor ^^(fc). It is propor- 
tional to the number Na{uj, k)dujdQk of photon states for a given polariza- 
tion a, frequency to, and direction k, since Na{uj, k)dujdVtk oc k^dkdVtk = 
nl^{k)uj'^du}dQk/c^- 

The equation of radiative transfer theory formally corresponds to a 
Boltzmann equation balancing all the changes in W^{R,T,t): 

= cY,I^B„,Ck,MriU9)WliR.T,t) + S^{R,T) . (10) 

The first term gives temporal variations of W^{R,T,t) due to e.g. time 
dependent light sources. The second term involves the divergence of the 
Poynting vector VgaW^{R,T,t). The correlation function W^{R,T,t) 
changes when there is a net flow of energy [t = 0) or correlation (t ^ 0) 
out of the volume element around R. The third and fourth terms describe 
losses and gains due to scattering. The scattering mean free path /^(fc) has 
been already introduced in Eq. (|^). Finally, S'^{R,T) indicates a source 
term. The transport equation ( p!OD is intuitively understandable. Com- 
pared to standard textbooks on multiple light scattering [27, 28] it contains 
two generalizations. First, it is valid for general anisotropic random media. 
Secondly, for t 7^ 0, it describes the transport of electric field correlations 
[11, 36]. The transport equation can be derived from first principles, i. e. 
starting from Maxwell's theory. It follows in a straightforward way from 
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the Bethe-Salpeter equation for the averaged "two-particle" Green func- 
tion [19, 20, 22, 29]. Its vahdity is restricted to length and time scales 
much longer than the wavelength and the time period of light. 

A. Diffusion Approximation 

To derive the diffusion approximation from Eq. (p!0[) we set t = and ne- 
glect any source terms. We first study the equilibrium solution of Eq. (]T0|), 
where W^{R, T) equals a constant Wq. As a result formula (0) for the scat- 
tering mean free path la{k) is reproduced. The energy density n^^{k)WQ 
still depends on the direction k and polarization a of light, due to the 
equipartition of the light energy on all available photon states. So far, 
this still needs an experimental confirmation. The diffusion approxima- 
tion follows when the k dependence of W^{R,T) deviates only slightly 
from the equilibrium angular distribution [27, 28]. Let us therefore look 
at an angular expansion of W^{R,T): 

W^{R,T) = ^Wo{R,T) + ^^nJ-J''{R,T) 

oTT 47r cri'^ 

+ Y: WmiR.TWJk) , (11) 

l>l,m 

where we introduced the total energy density 

^WoiR, T) = E / dQknl{k)W;^{R, T) (12) 

as Wq times an angular average over the cubes of both refractive indices: 

^= (;^ + ;^)/2 = y'dfifc[n?(fc) + r2^(fc)]/87r = (£y'£||+£5./')/2 . (13) 

Our goal is to establish a diffusion equation for Wo{R, T). We also defined 
the total energy density current J"'{R,T) of light with polarization a, 



J''{R,T) = I dn,nlvg^W^{R,T) 



(14) 
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A few comments are necessary to understand Eq. (fiA]) . In isotropic sys- 
tems, the second term on the right-hand side just corresponds to an ex- 
pansion into the components of = (sin'(9cos(/9, sin-i^sinyj, cos"!?), which 
basically stand for / = 1 spherical harmonics. Here, for uniaxial systems, 
it is useful to choose basis functions with nl^{k) as a weight function: 

J dcOsM^nl{k)Yi':^{k)Yi7^,{k)=^6u'Srnrn' ■ (15) 

We will see below that this choice establishes an approximation scheme 
for the diffusion constants of light. For ordinary waves we still obtain the 
conventional spherical harmonics. For extraordinary waves we introduce 
a new coordinate C = ni{k) cos{}/n2, which is equivalent to Uq - k = cos{} 
since it also ranges from —1 to 1. With this choice the weight function 
becomes a constant: 

J dcos^nl{k)... = ^J dC... (16) 

Hence, the basis functions Yil^{k) for a = 1 simply follow from spherical 
harmonics when cos is replaced by C. With the abbreviation C = cos 
for ordinary light, the generalized spherical harmonics are the same for 
a = 1 and 2. For Z = 1 in real representation, they read: 





In expansion ([Tl|) we have already used them explicitly in the second term. 
We also note that ni{k) sm.'d/n2 = ^e\\/e±\/l — C^. 

We are ready to extract the diffusion approximation from the transport 
equation (|10[). Multiplying Eq. ([T0|) by and summing over all directions 
of k and the two polarizations leads to the continuity equation for the 
energy density 

^-^Wo + V ■ J = . (18) 

The vector J = J^+J"^ denotes the total energy density current. A second 
equation, Fick's law, relates J to the gradient of the energy density: 
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J = -DVnm^ . (19) 

It will be derived below. We introduced the diffusion tensor 

D = D^1 + {D\\- D^)nQ®no (20) 

with its two independent light diffusion constants D\\ and respectively, 
parallel and perpendicular to the director tiq. Eliminating the current J 
finally gives the diffusion equation for Wq{R,T), 

[§f ~ " ^^^^) ^0^^'^) = ° ' ^21) 

where V = (V i , Vii). 



To obtain Pick's law (|T9D, we first look at components parallel to the 
director Uq. We project the transport equation (|^) on Y{^(k) oc Ck and 
arrive at a set of equations which couple the energy density Wq and the 
components J|| and J||: 



(47r)3 c ^ ( l \ J Bl \ ( J|| 

^12 ^22 



The quantitites BI^q are extended matrix elements of Bap{k, q): 



Bli = i I [{Cl-CuC,)B,,{krq) + cf-^B,2{krq)] 

Jk J q £\\ 



BI2 = L I [{Cl - C,C,)B,,ik, q) + cf-^B,,ik, q)] (23) 

Jk J q S\\ 



B\2 = - U CuC,B^2{krq) 

Jk J q 



where we used the abbreviation = / dCkdipk- In Eq. ( ^21) we neglected 
terms proportional to dJ^/dT. We also assumed that further contribu- 
tions containing Wiq with I > 1 are small and discuss this approximation 
below. Solving for the two currents J|| and J| leads to the parallel com- 
ponent of Eq. (|I^): 
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J|l = Jl + = -D||V||nWo (24) 



with the diffusion constant 



(477)3 ^ bI + bI^-2BI 



'12J 



The perpendicular component of Eq. ( [191) follows in an analogous 
way after projecting the transport equation (p!0|) on [l^"(fc) + Y{^_-^^{ik)] oc 



J± = Ji + Jl = -D^V±n^Wo (26) 



with 



^ 18 nl^ BhBi,~iBky • ^ ^ 

The matrix elements B^^ are defined as in Eqs. ([23| ) but with C replaced 

by Vl — C*^ COS(y9. 

The expressions ( p5|) and (|27|) for the diffusion constants are approxi- 
mate formulas since in Eq. (^2D we have neglected terms containing 
with / > 1. In ref. [20] we showed that corrections of D\\ and re- 
sulting from 1=3 spherical harmonics are essentially 1% or smaller. One 
could ask for the reason. The additional coefficients Wiq in Eq. ( p2D with 
the non-diagonal matrix elements {10\Bai3{k, q,0)\lO) as prefactors have 
to be determined by projecting the transport equation on higher gener- 
alized spherical harmonics (Z > 1). Because of the choice of our basis 
functions the resulting equations only couple to jj^. The coefficient 
Wq does not appear. That means Wiq (/ > 1) depends in a complicated 
way but directly on . As a result the matrix elements B^lj^ are renormal- 
ized. However, the renormalization is small if the non-diagonal elements 
{10\Bai3{k, g,0)|/0) are small. For symmetry reasons it is clear that D\\ 
and D± involve, respectively, m = or |m| = 1 spherical harmonics only. 
In isotropic systems higher spherical harmonics do not contribute at all, 
since {10\B{k ■ g,0)|/0) oc 5ii, and D oc {1 — cosi^^)"^ follows. 
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An extensive discussion of the diffusion constants and tfieir depen- 
dence on Prank constants Ki, dielectric anisotropy Ae, and applied mag- 
netic field H is given in refs. [20, 21] and by Tiggelen et al. [22]. We 
just summarize the important results. For a typical material, 5CB, we 
find D\\ = 1.43 X lO^cm^/s and D± = 0.98 x lO^cm^/s with a ratio of 
D\\/ D± = 1.45 in excellent agreement with experiment and numerical sim- 
ulations [21, 25]. The values were calculated in the limit of — > 0, which 
demonstrates that the divergence of the structure factor for — > does 
not affect the diffusion constants. If one introduces transport mean free 
paths via — 3n2D\\/±/c, one arrives at l\\ — 2.2mm and l± — 1.5mm 
for 5CB. We stress that it is not obvious how to define transport mean 
free paths in anisotropic turbid media. The calculated values give an ori- 
entation only. In the case of Ae = and the one-constant approximation 
Ki = K there is still a remaining anisotropy of D\\/ D±_ = 1.06 due to the 
inherent anisotropy in the nematic structure factor. For Ae < there 
exists a point with = D± beyond which the ratio D\\/D_[_ is smaller 
than 1. This behavior should be observable in discotic nematics. 

B. Dynamic Absorption 

Now, we set t 7^ and look at the transport of electric field correlations. 
In doing so, we restrict ourselves to times t much smaller than typical 
director relaxation times r = 7/ {Kq^). For wave numbers qs of light they 
cover the range r = 10 — 100 jis. We rewrite the dynamical structure 
factor, 

Ba0{k, q, t) = Bap{k, q, 0) + [Ba^ik, q, t) - B^p{k, q, 0)] , (28) 
where the second part on the right-hand side assumes the form 

B^pik, q, t) - B^fsik q, 0) = -{Aefk^T^ ^^^^^M t (29) 

5=1 Vsyqs) 

after an expansion of the exponential factor in B^pik, q,t). All the coeffi- 
cients in an angular expansion of W^{R, T, t) now carry the relative time t 
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as a further argument. The important quantity is the total autocorrelation 
function 



71' 



Wo{R,T,t) = J2 J dQknl{k)W^{R,T,t) 



(30) 



In repeating the derivation of the continuity equation (|T8|) we have to add 
a dynamic absorption term fi{t)n^Wo{R,T,t): 



nmo{R, T, t) + V ■ J{R, T,t) = , 



(31) 



which means that temporal correlations are not conserved quantities be- 
cause they decay to zero. The dynamic absorption coefficient fi{t) follows 
generally from an angular average over all dynamic modes of a system: 



(32) 

Further terms in Eq. (|3ll) containing coefficients are small, due to the 
assumption that Baf3{k, q,0) ^ Ba/sik, q,0) — Bap{k, q,t), and can be 
neglected. The same applies to Fick's law in Eq. (|19D, so that we arrive 
at the diffusion equation with absorption which we already formulated in 
the introduction. 

Using Eq. (^^ for director modes we calculate 

fi{t) = fxot with Ho = — ^ , — - , (33) 

97r y/6± 7 

where the numerical factor /2 stands for a dimensionless angular average 
over the geometrical factor N{a,P,S) and the viscosities r]s{qs)/l- The 
factor Jl is of the order of 1 in thermotropic nematics but always larger 
than 1 [20]. Note that /i(t) only depends on viscosities and not at all on 
elastic properties. They cancel because they determine both static light 
scattering and hydrodynamics of the director modes. For 5CB, 7 = 0.81 P 
which agrees very well with the experimentally determined value of 7//! = 
0.60 ± 0.20 P [21, 25]. In experiments one fits the autocorrelation function 
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with the help of a numerical or exact solution of the diffusion equation 
under the appropriate boundary conditions [7]. 

Down to the experimental resolution of 4 x 10~^ s no deviation of the 
director dynamics from the Leslie- Erickson theory was observed. It should 
show up in a different temporal power law of n{t). It would be interesting 
to study systems with higher viscosities like polymer liquid crystals and 
to look for such a deviation. The Brownian motion of colloidal particles 
in a fluid e.g. clearly does not show a simple diffusive behavior on very 
short time scales [7, 37]. 
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